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Abstract. We investigate the conditions for a bounce to occur in Friedmann - 
Robertson - Walker cosmologies for the class of fourth order gravity theories. The 
general bounce criterion is determined and constraints on the parameters of three 
specific models are given in order to obtain bounces solutions. Furthermore, unlike 
the case of General Relativity a bounce appears to be possible in open and flat 
cosmologies. 



1. Introduction 

The idea of a "bouncing" universe in Friedmann- Robertson -Walker (FRW) 
cosmologies has been examined many times since the 1930's when Tolman [I] proposed 
that closed (k = +1) FRW universes might re -expand after collapsing to a high 
density state in the future J. This subject has remained popular throughout the 
history of cosmology |2j and is currently a very active area of research, motivated 
primarily by recent developments in M-theory, braneworlds and quantum gravity |3j . 
In particular, the ekprotic model of the universe 0] is one such realisation of a cyclic 
cosmology inspired by M-theory while in loop quantum gravity the semi - classical 
Friedmann equations have correction terms that produce a bounce [H] • 

In classical cosmology such bounces are not possible in FRW models if the active 
gravitational mass is positive: that is, if p + 3p > 0. This is a direct consequence of the 
Raychaudhuri equation which is the fundamental equation of gravitational attraction 
H3EIIH]- On the other hand quantum fields and indeed classical scalar fields c/> can 
violate this condition Hence, such fields can in principle allow bounce behaviour 
in FRW models, however such bounces are difficult to produce in universes that grow 
large enough to be realistic: typically the probability of a bounce is of the order of the 
ratio of the minimum to maximum expansion size (scale factor) ITT] . Furthermore 
very anomalous physical behaviour can occur if these classical fields violate the reality 
condition <\? > 0, which is equivalent to requiring that the inertial mass density is 
positive: p + p > 0. It is also possible that ghost fields, i.e., fields that have a negative 
energy density are capable of producing a classical bounce ^3 • 

It is not only FRW models that are of interest. Smolin's idea of collapse to 
a black hole state resulting in re - expansion into a new expanding universe region 
[I'M ITlj suggests that the geometry of the universe at a bounce might be very 

J It was already known that this was not possible for k = and k = — 1 models. 
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different than the completely isotropic and spatially homogeneous FRW spacetimes; 
indeed, because some spatially homogeneous modes are unstable > a more general 
geometry might be expected. For example, in the scenario of black hole collapse and 
subsequent re - expansion, we might expect the geometry at the bounce to be that of 
a Kantowski - Sachs model E| > because this model has the same symmetries as 
the spatially homogeneous interior region of the extended (vacuum) Kruskal solution, 
that represents the late stage of evolution of an isotropic black hole when the matter 
can be neglected. 

Because of the growing interest in bouncing and cyclic cosmologies it is worth 
exploring the full phase - space of possibilities both in General Relativity and other 
models of gravity. Indeed classical bounces in Kantowski - Sachs have recently been 
discussed in ^B] within the context of standard General Relativity. 

In this paper we examine the conditions for a classical bounce to occur in models 
of gravity that have an "effective" fourth order action given by 



Our purpose is to use the 1+3 formalism in order to deduce the condition for which 
fourth order gravity admits a bounce. Such conditions will be given in terms of the 
parameter(s) for three specific forms of f(R): R n ,cxp(\R). 

In the last few years higher order theories have been proposed as theoretical 
models for solving the problem of cosmological acceleration ^3 EH]- The study of the 
bounce conditions for such theories allows us to determine if in fourth order gravity 
there is a connection between a bounce and the cosmological acceleration phenomenon. 
If this is the case, we would be able to use the bounce conditions to put constraints 
on cosmic acceleration and vice versa. 

This paper is divided up in the following way: in section 2 we give a brief overview 
of the fourth order theories of gravity f(R); in section 3 we derive the bounce equations 
in the case of FRW cosmologies for a generic f(R) theory; in sections 4-6 we specialise 
the discussion to the three models listed above. Finally in section 7 we present our 
conclusions. 

Unless otherwise specified, natural units (ft = c = ks = 8ttG = 1) will be used 
throughout the paper, Latin indices run from to 3. The symbol V represents the 
usual covariant derivative and d corresponds to partial differentiation. We use the 
(—,+,+,+) signature and the Riemann tensor is defined by 




(1) 




The Ricci tensor is obtained by contracting the first and the third indices 




Finally the Hubert-Einstein action in presence of matter is defined by 




(•5) 
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2. Preliminaries 

By varying the action Q, we obtain the fourth order field equations 

f\R)R ab - l -f{R)g ab = f'(R) cd {9ca9db - 9cd9ab) + , (6) 

where the prime denotes a derivative with respect to R. It is easy to check that 
standard Einstein vacuum equations are immediately recovered if f(R) = R. 

In the 1+3 covariant formalism |22], after choosing a frame comoving with the 
average matter velocity vector u a , the Raychaudhuri equation for a fourth order 
gravitational Lagrangian can be written as 

9 - V a u a = - 1 9 2 + (u a ii a ) - 2(a 2 — us 2 ) — -^-p (n + 3p) 

- (/ - Rf + 3/' + 9/' - A/') , (7) 
and the Gauss- Codazzi equation (see equation (55) in [22]) is 

e 2 + 3^9 = ^-^ + 3a 2 -i[(/-i?/') + A/'] , (8) 

where / = f(R), f'{R) = df /dR, R is the 3-Ricci scalar, A is the Laplacian operator 
and a "dot" corresponds to the covariant derivative along u a . Note that in this 
equation, as well as those that follow, we consider the Ricci scalar R as an independent 
field. Such a position was first proposed in canonical quantisation of higher order 
gravitational theories [2J and allows us to write the 1+3 equations as a system of 
second order differential equations at the price of adding the constraint 

i? = 2^e + ^e 2 + a 2 -2^ 2 + ^ . (9) 

If we now consider the trace of (J5J) 

3/' + 36/' - 3A/' = f'R - 2f + y. - 3p , (10) 
we can simplify the Raychaudhuri equation to give 



e - w a u a = - ± e 2 + (u a u a ) - 2(a 2 - u?) - jj - -L + jj a/' - e/' 



(ii) 



Equations (JSJ) and the l|ll|l are the key equations we need to determine the bounce 
conditions for this class of cosmological models. 



3. Bounce Equations 



In what follows, we define the occurrence of a bounce at time t = by the conditions 

9(t b ) = , 9(4) > . (12) 

In FRW models §, such conditions can be written as 

S(t b ) = , S(t b ) > . (13) 

§ In anisotropic models which have more that one scale factor Xi, i = 1,2,3, the 1121 has to be 
understood as characterising a bounce in the average scale factor of the universe S = \/ X1X2X3. 
However one can also consider a more generic situation where a bounce could occur in any of the 
scale factors X{. We can make this precise by defining the expansion parameters X4 = Xi/Xi, so a 
bounce in Xi will occur at time t = tj, iff Xi(ti) = and Xi(tb) > 0. It is clear that although it may 
be possible to have a bounce in one of the scale factors but not the other, this does not lead to a new 
expanding universe region. We therefore require that a bounce occurs in all Xi's, even though they 
may in general occur at different times. 
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Let us now apply the bounce condition to equations © and (|llfl . In the case of FRW 
models (where a — 10 — ii a = 0) at t = t b we obtain 

©* = -^ + A, R b = 2^ + ^(f b -R b fi) , (14) 

Jb Z h Jb Jb 

and 

K = 2( •■ — | - (I'.) 




where we have used the suffix b to indicate quantities evaluated at the bounce time 
t b . It is clear from the first of (|14|l that the presence of the non-minimal coupling 
affects the sign of Q b in two ways: through the non-minimal coupling of higher order 
terms with matter and the non-linear contribution to the Lagrangian. These two 
corrections are related to each other, but in principle a bounce can be induced by just 
one of them. The second equation leads to a surprising result: in higher order gravity 
a bounce does not necessarily imply positive spatial curvature. This does not happen 
in standard GR, which can clearly be seen when f(R) = R. 

Another interesting aspect of the system above, is that the second equation in 
l|14f> is not independent of the first one. This can be easily seen by substituting (|15fl 
and then solving for 0. For this reason, we will consider only the first equation of 
Ill4|l in our calculations: 

° k = -f + 4' (16) 

Substituting for the expansion in terms of the scale factor S and using the standard 
result for the 3 -curvature for FRW models: 

Q = 3S/S, R = 6K/S 2 , (17) 
equations I|15ll6fl can be written as 

n V-b . fb D „ ( S b K \ 

Qb = ^ = -Ji + M' R = 6 {s-b + sl) ' (18) 

Equations i|18fl are the defining equations for a bounce in FRW cosmologies and 
generalise the well known results for General Relativity to f(R) models. From now 
on, for sake of simplicity, we will continue using the variables (0b,i?b), noting that 
the sign of (0b, Rb) correspond to the sign of (S b , K). 

It is important to realise at this point, that since equations Ijl8|l are highly non- 
linear, we cannot find an exact solution for them. However, since our purpose is simply 
to investigate the possibility that a bounce may occur and not specific bounce models 
in fourth order gravity, we can simply limit ourselves to checking the consistency of 
the system ()18|l subject to the bounce conditions Ijl2|l . 

In the remaining sections of this paper we consider three specific forms of /(i?): 
R n , R + aR m and exp(Ai?). 

4. f(R) = R n 

In R n - gravity, the function / is specified by a generic power of the Ricci scalar: 
f(R) = R n |191 1201 \2'i\ . In our analysis, the trivial cases n = and n = 1 are 
neglected and n may either be a relative or a rational number. 
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Table 1. Sign of the 3-Ricci scalar for different type of values of n in R n gravity 
bounces. The numbers r, q are chosen to be relative (r, q S Z). 



n n > 1 n < 1 









n E Z odd 


R b >0 


i?b^0 


7i E Z even 


R b > if R > 


i?b ^ if i? > 




#b < if R < 


R b < if i? < 


2q 


R b > 




2r 

71 2q+l 


.Rb > if i? > 


7?b ^ if R > 




R b < if i? < 


R b < if i? < 


n - 2r+1 

" 2q+l 


i?b > o 


i?b^0 



In this case the bounce equation (|16(l is 
which can be written as 

<"-"*'— ^r + T ■ fl ^ 2 ( e ' + l)^ < 20 » 

In a previous paper, focused on the dynamics of R n gravity |23j . it was shown that 
the sign of the Ricci scalar remains unchanged once the initial conditions are fixed. 
This means that, if at the bounce we find the Ricci scalar to be positive, the sign of 
R is positive for the entire cosmological history. Hence, we can consider the sign of R 
as an "initial condition" . 

The form of (|20|l suggests the existence of two different behaviours depending on 
whether n is bigger or smaller than 1 (see Table If n > 1, the combination of (|19fl 
and HI 5(1 reveals that once the sign of R and the set of values of n is fixed, the sign 
of the 3-Ricci scalar is uniquely determined. In particular, a closed bounce occurs for 
most of the values of n and R with the only exception of n even with R < and n 
rational with odd denominator and even numerator with R < 0. The same cannot be 
said for the case n < 1. In fact, with the only exception of n even with R < and 
n rational with odd denominator and even numerator with R < 0, for which have an 
open bounce, the equation 1(20(1 admits both an open and a closed bounce. 

These results are very interesting if combined with those found in the dynamical 
analysis of [22|. Here it was shown that a cosmological history exists which contains 
an "almost-Friedmann" phase (AFP) followed by an accelerating phase (ACP). Such 
a history is possible only for 

n t lo / « < 0.37 if neAfeven, / 01 x 

W = °' 1/3 1 1.37<»<2 if neN odd , W 

( n < 0.37 if n € N eV en , / 00 x 
W = 1 1 L5<n<2 if neM odd . (22) 

where Af even is the set of even integers or rational numbers with an even numerator, 
and Af odd is the set of odd integers or rationals with odd numerator. 

These values of n can be traced back to our variables giving i?<0if0<n<l 
and R > if n < and n > 1. 
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Comparing this with our results (beware of the change in the metric signature!) 
and bearing in mind that the sign of R is fixed for a given cosmological history, we 
conclude that if n > 1 and n £ Af dd, the presence of AFP— >ACP orbits is compatible 
with closed bounces for cosmic histories with R > 0. The same happens with Af eV en 
and n < with i?>0or0<n<l with R < 0, where in the second case n can be 
only rational. 

Furthermore, the analysis of i?"-gravity with the supernovae type la given in 
|24| shows us that n must approximately lie in the range 1.366 < n < 1.376. This 
is within the limits given in equation 121fl and therefore the possibility of a closed 
bounce occurring in such cosmological histories is not excluded. 



5. f(R) = R + aR m 

This model has been one of the most popular fourth order gravitational theories, and 
is still among those most studied. There has been particular emphasis on the sub-case 
m = 2, because these corrections have been shown to stabilise the divergence structure 
of gravity, allowing it to be renormalised (at least at first loop) ■ In our study the 
parameter a will be taken to be real, but m can either be relative or rational. 
The general bounce equation l(TH|) for f(R) = R + aR rn reduces to 

■ R b + aR' b n -2^ b 

° b = 1 + maRr 1 ■ (M) 
Using the (|15|l the same equation can be also given in terms of the 3-Ricci scalar 
- a(m-l)iy + 2 M6 

R " = 1 + maR™- 1 ■ (M) 

For our purposes both these versions of the bounce equation are useful. In fact, from 
the first one we can obtain the values of the Ricci scalar for which a bounce is possible 
(i.e. Ob > 0) and then use them in l|24|) to obtain the sign of the spatial curvature. 
This is possible because the Rb term contains both and Rb and the two equations 
are non-linear in Rb, so that the system (|23I24[1 may be considered as a parametric 
version of the full relation Q b (R b ). Unfortunately, the complexity of the (|23|l and (|24|l 
makes it impossible to find general exact results. However, since most of the important 
features of the RHS of (|23|l and (|24|l (number and sign of the solutions, etc.) depends 
only on the nature of m, we can still give some general results. 

A comparative analysis of the two equations above leads to the results listed 
in Tables 2-3, where the quantities i?*, Ri, Rf represent the values of R for which 
Ob = 0, Rb = 0, and (1 + mai?™ -1 ) — > respectively. We can see that closed bounces 
are allowed for every integer value of m (often together with open bounces) . For m 
rational, closed bounces are not allowed in general for < m < 1. For m rational 
with even denominator we have no closed bounce for (m > 1, a < 0) and no bounce 
at all for negative m and a. 



6. Case /(i?) = exp (XR) 

As a third example, we consider theories with a Lagrangians which may be expressed 
as an exponential of the Ricci scalar. This type of Lagrangian is interesting because it 
contains in some sense the previous two cases due to the fact that the exponential can 
be developed in powers of the Ricci scalar. In other words, the study of an exponential 
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Table 2. Sign of the 3-Ricci scalar for m integer. The quantities R* , Ri, Rf 
represent the values of R for which Q b = 0, R b = 0, and (1 + mai?™ -1 ) — > 
respectively. The index i, when present, indicates the presence of multiple roots. 
The root with lower numerical value have lower index i . Even if the last two 
quantities can be calculate exactly in terms of the values of the parameters, the 
first one requires numerical calculation. 



(a, m) Q b (R b )>0 R b 



to 6 Z even 


to > 0, a > 


R$<Rb< R Q 
R b > i?J 


.Rb < 
7R b > 




TO < (T < 


7?° < Rh < 

-Tlfc > U 


R b < 

E? \ A 




to > 0, a < v 


it u < R b < R\ < 


R b < 
Kb > 




. ^ A „ ^ A 

to < U, a > U 


D* / > ^ A 
< /tfc < U 

i?2 < -Rb < i?i 
Ri<R b < R Q 
H b > rt 3 > U 


D ^ A 
-Rb < U 

vr 6 > o 

Kb < 
Hb > U 


to 6 Z odd 


to > 0, a > 


R b > i?* 


Kb > 




to < 0, a < 


R* < R b < 
i? 6 > i?* > 


R b < 

k 6 > o 




to > 0, a < 


Ri < i?b < i?® 
R^ < i?h < .R 


i? 6 < o 
k 6 < o 




to < 0, a > 


R b > i? 


k 6 > o 


V 

TO = ±~ 

2q 


to > 1 , a > 


-Rb > R 


K b > 


p,qe Z 


0<TO<l,a>0 


R b > it 


K 6 > 


p odd 


to < 0, a < 


no bounce 






to > 1 , a < 


VR 


K 6 < 




0<TO<l,a<0 


0<R b < R Q 


R b > 




to < 0, a > 


Ri < R b < -R 
R< R b < R Q 

R b > i?2 


k 6 > o 

K 6 < 
K b > 



Lagrangian gives us the chance to investigate, in a relatively easy way, what happens 
if we consider a Lagrangian made up of a combination of different powers of the Ricci 
scalar. In our analysis, the parameter A is taken to be an arbitrary real number. 
The general bounce equation J3J for an exponential Lagrangian becomes 

® b = 2A exp (~ Ai?b ) [ ex P( A ^) ~ 2 W>] • (25) 

This equation tells us that in order to have a bounce (Q b > 0), one of the two conditions 
below have to be satisfied: 

A>0, f A<0, 

R b > i R b < ^ . (26) 
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Table 3. Sign of the 3— Ricci scalar for m rational. The quantities R*, Ri, Rf 
represent the values of R for which b = 0, Ri, = 0, and (1 + mctR™" 1 ) — » 
respectively. The index i, when present, indicates the existence of multiple roots. 
The root with lower numerical value have lower index i . Even if the last two 
quantities can be calculate exactly in terms of the values of the parameters, the 
first one requires numerical calculation. 



(a, to) 



® b (Rb) > 



R,, 



2q+l 

r,q£ Z 



m > 1, a > 

< to < 1, a > 

to < 0, a < 

to > 1 , a < 

0<TO<l,a<0 
to < 0, a > 



Ri < R b < 
R* 2 < R b < 
R Q < R b < 

Rb > R 
R < Rb < RI 

R b > R* 2 
R 3 < Rb < R 
Rb > R 
R b >0 
Rl<R b <0 
R Q <R b < R-2 
R 2 <R b < R 3 
Rb > Rt 



R b <0 
R b >0 
R b <0 
R b >0 
R b <0 
R b >0 
R b <0 
R b >0 
R b < 
R b <0 
R b >0 
R b <0 
R b >0 



TO - 2r+1 
m ~ 2q+l 


to > 1, a > 


Rb > Rt 


R b >0 


r,qe Z 


0<to<1,q;>0 


R b > Rt 


R b >0 




to < 0, a < 


R* 1 <R b <0 

Rb > i?2 


R b <0 
R b >0 




to > 1, a < 


i?t < i? fc < 
< i? fc < Ri 
Rb > i?i 


R b <0 
R b <0 
R b >0 




0<TO<l,a<0 


i?b > R-t 


R b <0 




to < 0, a > 


Rf < R b < 
i?^ < i?b < i?i 
i?b > i?i 


R b <0 
R b < 
R b > 



Taking into account ©, we see that only in the first case a spatially closed bounce is 
possible. This result was expected. In fact, if we perform a Taylor expansion of the 
function exp (Ai?) we have 

exp(Ai?) = 1 + XR + X 2 R 2 + .... , (27) 

so that in situations where the curvature is small, we obtain 

exp(Ai?) ra A Q + R}j , (28) 

which is equivalent to the Hilbert Einstein action with a non-zero cosmological 
constant. Consequently, A > implies a positive cosmological term in the action, 
which has been shown many times to lead to a closed bounce. 
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7. Conclusion 

In this paper we have derived the general bounce conditions for Friedmann Robertson 
Walker models in fourth order gravity and applied them to three specific models: 
f(R) = R n , f(R) = R + aR m and f(R) = exp(Ai?). 

Our analysis shows that in i?™-gravity a closed bounce occurs for most values of 
n and R with the only exception being n even with R < and n rational with odd 
denominator and even numerator with R < 0. Comparing these results with the those 
obtained in |2.'{| , we obtained the values of n for which a cosmological history consisting 
of an almost-Friedmann phase followed by an accelerated expansion is compatible with 
those that give a closed bounce and found that they lie in the same interval compatible 
with observations of type la supernovae, given in |24j . 

The R + aR m models are more complicated due to the fact that there are two 
parameters involved. Closed bounces are allowed for every integer value of m (often 
together with open bounces), but this is not true for rational values of m. In the case 
m = 2, a > 0, our results are in agreement with the independent study done by Page 
PS] using the Gibbons-Hawking-Stewart canonical measure. 

For exp(Ai?) models we found that a closed bounce occurs if A > and the value 
of the Ricci scalar at the time of the bounce is a particular function of energy density 
at the time of the bounce. 

It is important to stress that the constraints found above are necessary but not 
sufficient conditions and therefore this paper only addresses the possibility of a bounce 
and not a specific bouncing cosmological model. Nothing in our discussion implies 
directly that the bounce occurs at an early stage of the Universes' history^. 

The most striking result of our study is that, contrary to what happens in 
standard General Relativity, a bounce is possible in cases in which the cosmology 
is not spatially closed. The idea of a bounce in an open or flat universe may appear 
difficult to visualise, but can be understood if we remember that the quantity 0& gives 
a measure of the deviation of matter worldlines. In this sense the bounce condition 
l|12|) simply means that there exists a phase in which the separation between the matter 
worldlines decreases to a minimum and then increases again. Since this phenomenon is 
independent of the spatial geometry of the spacetime, the bounce itself is independent 
of it. On the other hand the existence of a k < bounce indicates that the problem 
of the cosmological bounce in fourth order gravity cannot always be related to the 
problem of gravitational collapse in such models. Indeed, if we compare the Friedmann 
equations in General Relativity 



• • . •"/< % ■ (29) 
and in Fourth Order Gravity 

62 - 3 6 r 9m + wm im - Rfm ~ m ~ £ = ■ (30) 

we realize that, because of the non-linearity of l|3UI) many switches between contraction 
and expansion are in principle possible. In other words, the cosmological history may 
be characterised by a sequence of contraction and expansion phases, whose number 
and duration cannot be determined without an exact solution of the cosmological 

1 One could always argue that higher order corrections become important when the curvature is high, 
typical of the very early evolution of the Universe, but in this case we assume that our modification 
to gravity holds in any curvature regime. 
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equations. Such "wobbling" could in principle be detected by its effect on structure 
formation and CMB anisotropies, and in this way provide useful observational 
constraints on higher order gravity theories. 

How is it then possible to determine if a bounce "d la Tolman" , separating two 
more or less causally disconnected phases of the universe, is possible in fourth order 
gravity? A desirable constraint comes from the fact that we want such a bounce to be 
an absolute minimum of the scale factor and that at this minimum the energy density 
is higher than what it is during all the processes (recombination, etc) of classical 
cosmology, but lower than what it is when quantum gravity effects become important. 
Unfortunately, there seems to be no way to understand if the different bounce solutions 
found above satisfy these criteria without using an alternative methods or performing 
a direct analysis of the equations. Consequently, the existence of a cycling universe 
remains an open problem in fourth order gravity. 
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